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Abstract 

The semiclassical collapse of a sphere of quantized dust is studied. A Born- 
Oppenheimer decomposition is performed for the wave function of the system and the 
semiclassical limit is considered for the gravitational part. The method of adiabatic 
invariants for time dependent Hamiltonians is then employed to find (approximate) 
solutions to the quantum dust equations of motions. This allows us to obtain cor¬ 
rections to the adiabatic approximation of the dust states associated with the time 
evolution of the metric. The diverse non-adiabatic corrections are generally associated 
with particle (dust) creation and related fluctuations. The back-reaction due to the 
dominant contribution to particle creation is estimated and seen to slow-down the 
collapse. 
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1 Introduction 


The canonical quantization of general relativistic isotropic systems carried out in suit¬ 
ably chosen variables leads to the whole dynamics being determined by the Hamiltonian 
constraint of the Arnowitt, Deser and Misner (ADM) construction. Such an approach is 
particularly useful if one wishes to study the semiclassical regime of a system of self grav¬ 
itating matter and has been applied to the collapse of a sphere of classical dust with the 
associated scalar field being related to time [1]. In a previous paper [2], two of the present 
authors again applied the ADM formalism to a model describing the collapse of a sphere of 
homogeneous dust [3] leading to a black hole and examined how Hawking radiation arises 
within such an approach. In particular the Born-Oppenheimer (or adiabatic) approxima¬ 
tion for the coupled matter-gravity system was consistently implemented according to the 
original formulation of Refs. [4, 5] (for alternative views see also [6, 7, 8] and for their 
comparison see [9]) by assuming the gravitational degree of freedom evolved slowly with 
respect to the matter degree of freedom. In this note we shall study corrections to this 
approximation. 

The classical collapse of an isotropic homogeneous sphere of dust [3] has been studied 
extensively and is treated in numerous text books {e.g., [10]). A particularly interesting 
aspect of such a collapse is that the interior of the sphere is a three dimensional space of 
constant curvature whose radius depends on time (in the language of cosmological models 
it is a section of a Friedmann universe). Further, it has recently been shown that boundary 
effects are in general absent for classical fluids with a step-function discontinuity of the 
kind associated with such models [11] (see also [12] for a more specific treatment of the 
sphere of dust). The treatment of a scalar field in a cosmological context is known and for 
a particular regime one has that the homogeneous mode of the free scalar field is related 
to dust, that is a fluid with constant density and zero pressure [13, 14]. One may then 
naturally ask to what extent is this possible for the spherically symmetric collapse of a 
massive scalar field, that is what is the effect of the boundary of the sphere. Obviously if 
the radius of the sphere is infinite one must reproduce the cosmological models. Therefore, 
whatever the difference is between the two cases, it must vanish as the radius of the 
sphere tends to infinity. Since a free massive scalar field is localized within its Compton 
wavelength and just “feels” objects at such a distance, the boundary will only affect fields a 
Compton wavelength away and its net effect on the sphere of matter will be proportional 
to the ratio of the Compton wavelength of the scalar field to the radius of the sphere. 
The former ratio is expected to be small for sufficiently large spheres (thus 3-curvature is 
relatively small), a condition which has been previously noted and related to a breakdown 
of the adiabatic approximation for dust [2] or to the validity of the semiclassical (WKB) 
approximation [15]. 

Let us illustrate the above consideration by a simple model. In the comoving frame 
the dust particles are at rest and the adiabatic approximation implies that the radius of 
the sphere can be kept approximately fixed [2] on solving the matter equation of motion. 
Let us then consider a static spherically symmetric space-time metric 

-|- dr^ -|- dD^ , (1.1) 

which of course corresponds to a flat 3-space and also agrees with the Schwarzschild 
vacuum far from the event horizon. The action for a spherically symmetric scalar field 
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^(r, r) with inverse Compton wavelength ^ will be given by 


5$ = 47r J dr J dr - - 2V 


( 1 . 2 ) 


where ' = d/dr, ' = d/dr and we have introduced a potential V which is responsible for 
the spherical confinement of the scalar field and is zero inside the sphere (r < ro). The 
equation of motion for is given by 


-$' + ^^$ + 21 / 4 > = 0 . 


(1.3) 


Since we seek an analogy for dust, we want a solution that is homogeneous for r < ro and 
zero for r ^ ro (no matter outside the sphere). Such a solution is given by 


$ = (^(r) 


fir) 


(1.4) 


with (j) + (j) = 0, since our solution is stationary, and 


/ = 


r r < ro 

ro e~r (^-'’o) r > ro 


(1.5) 


in which we have made the choice oi f jr dimensionless (so that cj) has the same dimensions 
as <h) and any other dimensionless factor has been absorbed into (p. The corresponding 
potential is then for r > rg 


2 ’ 


( 1 . 6 ) 


plus a Dirac delta singularity at r = ro which ensures the continuity of the derivatives of 
/. One may second quantize 4>, by introducing creation (a^) and destruction (a) operators 
for the above solution, obtaining [m^ = h fi is the mass of a scalar quantum) 




V 2 

L -1 


/ 


(1.7) 


and subsequently evaluate quantities of physical interest for any state of dust quanta. One 
finds for the energy density [13, 14] 

p = i (($)2 + (|>')2 + ^2 |,2 ) 


// (a'l' a + 


rritp {d^ d + ^ 




l2' 


r < ro 


(l + ^ + j^) r>r„, 


( 1 , 8 ) 
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which corresponds to a constant density inside the sphere and an exponentially decreasing 
density with range ~ 1/// outside, as desired. Lastly one may evaluate the pressure P 
[13, 14] obtaining 




r < ro 


—m^ jj,^ {a^ a + ^ 


2 ^0 „-2 u(r—rn) /i i 1 i 1 A ^ 

3^e [^ + 2jrF + TJ^) r>ro 


(1.9) 


which is zero inside the sphere and its magnitude decreases exponentially outside. 

We may at this point, as in cosmological models, obtain an effective action for 4>{t) by 
substituting the particular spatial function /(r), Eq. (1.5), in the original Eq. (1.2) and 
integrating over the spatial coordinate r: 



Thus, unless one has tq ~ 1/^, the edge effects are negligible as we heuristically indicated 
before. Of course in the above simple model we have not considered all possible modes for 
the scalar field, but only the mode corresponding to dust and we have seen that for such 
a mode upon quantization one obtains, for any state, pressureless dust inside the sphere 
[ 2 ]. 

In section 2 we briefly review the effective action used to describe the collapsing sphere 
of dust along the lines suggested above (which will allow us to maintain the Robertson- 
Walker form of the metric) and the results obtained in Ref. [2], in order to prepare the 
ground for the study of the back-reaction on the geometry induced by quantum evolution 
(collapse) of the dust. This will require relaxing the adiabatic approximation and is the 
aim of the present note. 

In section 3 the method of adiabatic invariants [16, 17] is illustrated and is applied 
to solve the matter Schrodinger equation with a time dependent Hamiltonian obtained 
through the Born-Oppenheimer decomposition of the matter-gravity wave function and in 
the semiclassical limit for gravitation. After the identification of the parameter describing 
the adiabatic limit and the associated states, coherent states for the matter wave function 
having correct classical limits are constructed and used to determine the expectation value 
of the diverse quantities of physical interest. Since the collapse is adiabatic at some initial 
time, from the expressions obtained we estimate and compare the subsequent deviations 
from the adiabatic approximation and see the effect of back-reaction on gravitation. Lastly 
in section 4 our results are summarized and discussed. 

We use units for which c = 1, k = SttG, the Planck length is then £p = y/nh and the 
Planck mass is rup = h/ip. 
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2 Semiclassical collapse of a sphere of dust 


As mentioned in the Introduction, in order to describe the evolution of a collapsing sphere 
of homogeneous dust in vacuum it is convenient to consider a time-dependent scalar field 
(j) confined inside a spherical portion of a Robertson-Walker space-time (the interior of the 
sphere) with line element 


ds^ = 


-drf + —2 + dVP' 

1 — ep^ 


0< P< Po , 


( 2 . 1 ) 


where e = 0, ±1. The radius of the sphere then follows from the matching condition with 
the external Schwarzschild metric and is given by ro = pqK [2, 10]. The corresponding 
total effective action (gravity plus matter) is given by [18] 


S = 


1 

2 



(kk"^ -eK^ +K^ ((^2 


( 2 . 2 ) 


where dr = K dp is the proper time of an observer comoving with the dust. 
The Hamiltonian obtained from the action in Eq. (2.2) is 


H = - 





(2.3) 


where = —KK/k and (j). Canonical quantization leads to the Wheeler- 

DeWitt equation, 


1 

2 


K 


1 


--K - 




+ (j? K 


T(i^,,/.) = 0 , 


(2.4) 


where we have chosen a suitable operator ordering in the gravitational kinetic term [2]. 
One now expresses T in the factorized form 'k(A', 0) = K which, after 

multiplying on the l.h.s. of Eq. (2.4) by x* integrating over the matter degrees of 
freedom, leads to an equation for the gravitational part [2], 


1 

2 




52 


- - K" + 


1 


y dK"^ K j 


(x I x) 


(x 


iX2 




X/ 


Ip 


d 


ix\i^ 1- 


2(xlx)''"'5iX 


x)(x 


(xlx) J dK 


T^lx) 


, I 52 ~ kP ? 

-(x li^l x)V' = —^ 


2{x\x)^"'dK‘^ 


92 

dK"^ 


^ , 


where we have defined a scalar product 


(2.5) 


(X I X) = J d(px*i<P,K)xi<P,K) , 


( 2 . 6 ) 


and we have set 


V, = e-* rMK')dK'^ 


X = e^^r^^K')dK'^ 


(2.7) 
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If we now multiply Eq. (2.5) by x and subtract it from Eq. (2.4) we obtain [2] 
7.. if. J , ...2 1 / \ 




dK^ / dK^ 


dK^ \ dK^ 


dK / \ dK 


which is the equation for the matter (scalar field) wave function. 

On neglecting the r.h.s. of Eq. (2.5) one may introduce time [19, 4, 5] by writing a 
semiclassical (WKB) approximation for the wave function ip: 


dSeff 

dK 


e+i^eff 


( 2 . 10 ) 


where Seff is the effective action satisfying the Hamilton-Jacobi equation associated with 
the l.h.s. of Eq. (2.5), 


S,ff = - 


— [ dr ' 

2k J ( 


K,kl-{eK,-2K{H, 


( 2 . 11 ) 


and {Hm) is now evaluated for K = Kc which is where ip, Eq. (2.10), has support. One 
may then define a (proper) time variable. 


dip d 

'dk dk 


h dK 2 \ dK J dK^ 


_ iK~d I fdS^ffY^ d^S^ffjd 


Kh^ dr 2 \ dK 


dK^ ^ dK 


( 2 . 12 ) 


Further if the r.h.s. of Eq. (2.9) and the second term in the r.h.s. of Eq. (2.12) are 
negligible, one obtains the Schrodinger equation. 


. ^ dXs 1 I 2 . 2 7^ 3 i,2 


where we have scaled the dynamical phase 


Xs = Xs 


Xs = X exp 




(2.13) 


(2.14) 


and omitted ip while setting K = Kc, which is where the semiclassical gravitational wave- 
function has support. 

In Ref. [2] Eq. (2.13) was solved by making the adiabatic approximation Kc/Kc <C 
1 and evaluating {Hm) for such solutions. Eq. (2.11) then becomes the action of the 
Oppenheimer-Snyder model [3], 


= - ^ I dr [Kci kl - (e Kd - 2 RTq)] , 


(2.15) 
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with Kq constant and Kd is one of the classical trajectories 


r Kd = Kod^hM 

\ r = Koheir]) 


heiv) = { 


Tj — sin 7] 

e = +1 


rf /Q 

e = 0 

(2.16) 

sinh rj — rj 

e = -1 . 



The matching condition between the inner metric and the external Schwarzschild metric 
then gives for the mass parameter of the sphere 

M/k = Kqpo^/k = (H^) pq^ = Nm^po^ , (2.17) 

where N = {Hm is the (constant) number of dust particles. 

Let us emphasize that Eq. (2.17) introduces the size of the ball in the model. If, as 
suggested in the introduction (as well as in Refs. [2, 15]), boundary effects are negligible, 
the above is the only difference with respect to cosmological models where there is no 
exterior metric. Further the adiabatic approximation amounts to ( Hm ) = Hq pQ being 
constant, that is the number of dust particles is constant, although the frequency tOc = 
Hci { Hm ) of the corresponding Schrodinger state in the conformal time rj is not [2] (of 
course in the proper time the frequency = ( Hm ) is constant). 

In Ref. [2] it was subsequently verified that all the approximations were consistent for 
the collapse of the sphere up to its horizon radius = 2M the Compton wavelength 
of the scalar field is much smaller than the Schwarzschild radius of the sphere, 

2M>-, (2.18) 

which is precisely the condition prQ>2pM^l that one needs in order to neglect edge 
effects [15] (see Eq. (1.10)). 

In the following section we shall obtain solutions to Eq. (2.13) without making the 
adiabatic approximation, thus allowing for a change in the number of dust particles, and 
shall estimate the corresponding corrected classical trajectories Kc [20]. 


3 Adiabatic invariants and quantized matter 


A suitable method for the study of time dependent quantum systems is that of adiabatic 
invariants. In particular given a time dependent Hamiltonian Hm{t), a Hermitian operator 
I{t) is called an adiabatic invariant if it satisfies [16, 17] 


in^^ + [i{T),HM{r)] = b. 
ar 

(3.1) 

The general solutions to the Schrodinger equation, 


•t; dXsir) TT ( / N 

=Hm{t)xs{t) 

(3.2) 

are then given by 



(3.3) 


n 
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where | n,T)j is an eigenvector of /(r) with time-independent eigenvalue A„ and the {Cn} 
are complex coefficients. The phase (fn is given by 


Vn{T) = - I l{n,T' \hdr' +iHMiT')\n,T')jdT' , 


(3.4) 


/TO 


and is the sum of the geometrical phase whose associated connection is given in Eq. (2.8) 
and the dynamical phase displayed in Eq. (2.14). 

The Hamiltonian in Eq. (2.13) corresponds to a harmonic oscillator of fixed frequency 
fj, and variable mass K^.^. In this case it is useful to introduce the following linear (non- 
hermitian) invariant 


4(r) = 6(r) 


where the phase 0 is given by 


0 (r) = I 


dr' 


and the operator b by 


6(r) = 




0 ’ 


+ i {xTT^ — X Kc 


(3.5) 


(3.6) 


(3.7) 


The function x = x{t) is to be determined as a solution of the nonlinear equation 

I o • I 2 1 

X + ^ — X + ^l X = , 

Kc Kc x^ 

with suitable initial conditions. 

The system admits an invariant ground state (vacuum) | 0, defined by 

Ht) I 0,r)f, = 0 , 

and one can define an invariant basis of states H = {| n, r)^} with a Eock space 

(p)n 

Vn\ 

- e-®i^|0r) 

vn! 

where ipo may be replaced by —0/2 and 

b\n,T)f^ = ^/n\n-l,T)f, 

6't' I n,T)j = \/n + 11 n-b l,r)f, . 


(3.8) 


(3.9) 


(3.10) 


(3.11) 
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Thus, since [6, 6^] = 1, in the following we will refer to b and its Hermitian conjugate ¥ as 
the invariant annihilation and creation operators, and introduce the hermitian quadratic 
invariant 


h = n [ilh + -] =n b^b + 


1 


(3.12) 


with b the invariant number operator. 

The particle annihilation and creation operators a and al are defined by 


a(r) = 


I pK,- 

2h 


+ i 




pK.H ’ 


and the Hamiltonian operator is then 

Hm{t) = hp 




(3.13) 


(3.14) 


where iV = al a is the particle number operator which counts the number of quanta of the 
scalar field (j). One also has a corresponding vacuum | 0, r defined by a | 0, r )^ = 0 and 
a complete set of eigenstates A = {\ n,T )^}, 


such that 


\Jn\ 

a\ n,T)^ = Vn| n-X.r)^ 


at I n, r )^ = y/n + l | n + 1, r )^ , 


and [a, dt] = 1. 

The two different Fock space basis are related. Indeed one has 


(3.15) 


(3.16) 


from which 


where 





4> = \l\ x(b + bp 
i (5t_S) + xiF,3 (6 + 6t) 


r a = H*6 + H*6t 
[ at = H6t + H6 , 

.4M = iy;T7(x-^^-,:i) 

B(r) = l^(x + ^-,£) 
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(3.17) 


(3.18) 


(3.19) 



are the Bogoliubov coefficients. From the above one has an inverse relation, 


r b = Ba-A*a^ 

\ = B* o) - Aa . 


(3.20) 


The two basis A and B will coincide in the adiabatic limit for which the derivatives of 

o /o 

xKc are small. Let us suppose that for t = the two sets {| n, Tq )^} and {| n,T(, )f,} 
coincide. This is achieved if 


6(rJ = a(rj ^ h4(T-o) = i^M('ro) , 


(3.21) 


which leads to the following initial conditions for the function x 


x(rj = //iLc^(rJ 


- 1/2 


x(rj = 0 . 


(3.22) 


In particular, the second condition in Eq. (3.22) guarantees that the state of the system 
satisfies adiabaticity in the limit r ^ Tq. 

In order to study deviations from adiabaticity it is convenient to introduce a = 
X whereupon Eq. (3.8) becomes 




(T + 1 - 


3iLe 


\ 


fJ = 


<T + fj 



(3.23) 


As a hrst approximation [2], which we shall return to when we consider the effect of back- 
reaction, we set Kc which amounts to taking ( Hm ) as constant and equal to K^j k. 

This may be assumed to occur for the time tq for which the adiabatic approximation is 
valid (a and b coincide). One then obtains 


n = 


1 + 


3 e 5^ 


4 (d^he 


1/2 


1 + 


SeS^ 

8 {drjhe] 


(3.24) 


where <5 = In the above, 6 plays the role of an adiabaticity parameter which 

connects the departure from adiabaticity to the time as the collapse proceeds (r < tq). 
Indeed the collapse is adiabatic for 5 ^ 0, in agreement with Eq. (2.18). For the case 
e = 0, one has 0 = 1 and the exact solutions of the non-linear Eq. (3.23) (with Kc = Kd) 
can be obtained from the expressions displayed below by setting e = 0. For a general value 
of e (= 0, ±1) the solutions to Eq. (3.23) can be expanded in positive powers of 5“^ [21], 


OO 


cr = ^ (5^” Gn 

n=0 


(3.25) 


where Gn may be expressed in terms of O and its derivatives with respect to r. To the 
lowest order in (5^ one has 


G = + 0{5^) 


3e 

IQfJ^Kd^ ■ 


(3.26) 
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Further, given a particular solution a for a (and correspondingly for x), the general solution 
to the classical equation for cj) obtained from Eq. (2.2), 


i^c + 3 (/>£ + /i^ (/>c — 0 , 

J\r 


(3.27) 


is [22] 


(j)c = —p Kc ^ {D cos & + E sin 0) , (3.28) 

where D and E are constants and 0 is given in Eq. (3.6). We observe that this solution, 
obtained here through the adiabatic expansion of Ref. [21], actually coincides with the 
WKB-type solution of chapter 3.5 of Ref. [23]. In fact, the non-linear equation (3.103) 
in [23] can be recovered from our Eq. (3.23) by defining W = 1/cj^ and the vacuum 
corresponding to a as given above is the second adiabatic order vacuum in the terminology 
of Ref. [23]. 

The introduction of the eigenstates of adiabatic invariants allows one to consider the 
effect of particle production due to the variation of the metric on the evolution of matter. 
In order to see this, it is convenient to introduce coherent states in the b modes. 


lbs ® 


- ,-l«l72 


oo r. 

2/0 a'' 


E 




= e 


oo ri i n © 

_ ^-|a|2/2-i©/2 ^ a e 


n =0 


E 

n =0 


\/ri\ 


n,T)f, , (3.29) 


where a = u -|- f u is an arbitrary constant. We note that the modulus squared of the 
coefficient of j n, r in the above equation satisfies a Poisson distribution with a maximum 
at n = jaj^, average value and standard deviation jaj. If for some time Tq the adiabatic 
approximation is valid then the above Eq. (3.29) will also be a coherent state in the A 
basis. However, subsequently it will correspond to a squeezed state in the A basis, since 
the B basis is related to the A basis through a Bogolubov transformation Eq. (3.20). 

One may evaluate the expectation value of the diverse quantities of physical interest 
with respect to the above state. On defining (O = (a, r j O j a, r one obtains 


(j)c ^ = V2h jaj X cos(0 — P) , tan P = vju , 

in agreement with Eq. (3.28). Similarly for the momentum one obtains 


7r</,,c = 


a 


Kfp X cos(0 — P) — — sin(0 — P) 

X 


Eurther, 




( )bs = ( )bs = ^ 

both of which are independent of a and lead to 

{(Air C (>;/=5 /tErTTF . 


(3.30) 


(3.31) 


(3.32) 


(3.33) 
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Again this result is independent of a and we note that the uncertainty relation above is 
a minimum in the adiabatic approximation when the two Fock spaces A and B coincide. 
One also has 


i^c)bs — ^ \ + 


(3.34) 


{HM{T))bs = ^ 

^ 2 ^ ^ ^ ^ {^)bs (3-35) 


and, of course, in the adiabatic approximation, 

fJ'{Ic)bs-{HM{T))bs ■ (3.36) 

It is of a particular interest to examine the behaviour of the above quantities when 
one has (small) deviations from the adiabatic approximation. This may be achieved by 
considering corrections of 0{5‘^) to the adiabatic approximation and will allow us to see 
the effect of matter not following gravitation adiabatically and the corresponding back- 
reaction on gravitation. In particular one has 

2 fj 1 r 0x2' 

(t)r. ^ lal 1 - . „ . .o cos(0 - 0 ) , (3.37) 

H (Ko dr^ 


3 g 1 


where now 


Further, 


1 n 

Qc^- / 
d Jnci' 


f dC,dc^he 1 -h 
Jm{To) 


3 e 

4(9cM' 


/2,3 ( ^ 36(5^ . 


3 58"^ he 

2 (5„ heY 


cos(0 — /3) > , 


(3.37) 


(3.38) 


(3.39) 




((Avr<^)2 )bs - Ys ^0 1 ^ + 




3 e 5 

)3 [ S{d^he)\ 

(3.40) 

3(52 Tg 

^{dr^heYyl {dr^heY\] 

(3.41) 

\ 9 

_ 4 {d^he)\ • 

(3.42) 
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Also, the expectation value of the matter Hamiltonian is given by 

~ ^ + sin(0-/3) + ^(5^^^^^cos(0-/3) +cos^(0-/3)| 

/ 1 \ fi 

= ( |ap + - j + 3 h/U |ap(5 ^ ^ ^ sin(0 — (]) cos(0 — /?) 

= -^ + A{HM)bs^ (3-43) 

where only the first term K^jn on the r.h.s. of Eq. (3.43) survives in the adiabatic 
approximation and the second term is associated with particle production. 

In order to better understand our results, in particular Eq. (3.42) and Eq. (3.43), we 
may consider rj^ small and ;§> 5 corresponding to <C r <C tq, which implies that 
we are almost adiabatic (we only consider small deviations from adiabaticity). It is then 
straightforward to see that, in such a case, from Eq. (3.42), 

1 + 0 ($) . ( 3 ++) 

which of course shows that as i] becomes smaller (during the collapse), matter becomes 

less and less classical. On the other hand, from Eq. (3.43) one obtains 

, (3.45) 

again showing that, as rj becomes small, one has an increasing production of matter. 
Clearly such a production will induce a back-reaction on gravitation and Kc will change 
[24] from to + AK^, where AK^. is KoO{5/r]^). Since we are considering an 
exterior Schwarzschild metric and the time at which the horizon is crossed by the last 
shell of matter is given by 2M = po Kci'qn), the presence of the additional term AKc 
will lead to po Kciipn) — 2M — po Kq 0{5/pjj). This of course implies that the horizon is 
crossed at a later time (smaller value of p) and that the back-reaction (matter production) 
slows down the collapse. The effect of particle creation on a Friedmann-like collapse and 
the possible avoidance of the cosmological singularity due to quantum effects have been 
studied in Ref. [14]. 

In obtaining the above results we have neglected the r.h.s. of Eqs. (2.5)-(2.9) which 
are associated with fluctuations in the particle production due to the variation of the 
metric and the term 0{fi) in our introduction of time in Eq. (2.12). Let us estimate these 
effects and compare them with particle production. For this purpose, order of magnitude 
estimates of the r.h.s. are sufficient. They can be obtained by just evaluating the r.h.s. 
for a state | n = N such that N = \a\‘^, since it is for this value that the (modulus) 
coefficients of the expansion in Eq. (3.29) are peaked and corresponds to the amplitude of 
the classical oscillations (see Eq. (3.30)). For the r.h.s. of Eq. (2.5) one obtains 

nh^{N\^{l-\N){N\)^\N) = (^\Hm\L){L\Hm\N) 

^ L^N 

cs K,iA{HM)bsO(^^) , (3.46) 
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whereas for the modulus of the r.h.s. of Eq. (2.9) one has 




dK^ 




\N) 


E 

L^N 


l: 


n{N - L) 


{L 


d‘^Ic 

dK^ 


N)+ HL-P){L\Hm\P){P\Hm\N) 

P^N,L 




(3.47) 


Further, the corrections of 0{h) to Eq. (2.12) lead to an additional term for the 
Hamilton-Jacobi equation associated with Eq. (2.5), 


K 


3 irl 


ttk ^ Kttk 


4:K‘^7:t 2 K‘^'Kk 2K^ttk 


1 


a 


KaA{HM),p(^ 


rj- 


(3.48) 


whereas the modulus of the corresponding contribution to the r.h.s. of Eq. (2.9) is (omit¬ 
ting il>) 


kT%^ 


PJL. fA_(jv|A|jv 

2K7rK \dK ^ 'dK' ^ 


N] 


c^K,iA{HM)bsO 


(3.49) 


which is again of higher order with respect to KciA{ Hm )hs- 

The above results Eqs. (3.46)-(3.47) are associated with particle production due to 
the metric variations. In particular we see from Eq. (3.44) that, as the collapse proceeds 
the evolution of matter becomes less and less classical and one has the production of 
particles (Eq. (3.45)) which, as a consequence, slows down the collapse (back-reaction). 
Further one has lower order contributions, Eqs. (3.46)-(3.47), which are associated with 
fluctuations in particle productions and affect both the evolution of gravitation and matter 
and corrections to the classical gravitational motion (see Eq. (3.48)). 


4 Conclusions 

The matter-gravity system lends itself to a study analogous to that employed in molecular- 
dynamics where one also has two mass (or time) scales. We previously applied this analogy 
to the collapse of a sphere of dust in the adiabatic approximation, with the matter equation 
of motion being solved while freezing the gravitational degrees of freedom. In this note 
we have generalized such a treatment, to allow for the time dependence induced in the 
matter Hamiltonian by the time variation of the gravitational (metric) degree of freedom. 

Time has been introduced by considering the semiclassical (WKB) limit for gravitation 
while the matter equation of motion was solved by neglecting fluctuations and employing 
the method of adiabatic invariants. Such invariants, in contrast with the Hamiltonian, 
have time-independent eigenvalues. In particular, in our case, which corresponds to a har¬ 
monic oscillator with time dependent mass, on freezing the gravitational degree of freedom 
one reproduces the adiabatic results [2]. The adiabatic invariants are then immediately 
related to the usual harmonic oscillator (particle) creation and destruction operators and 
its Hamiltonian. 
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In general the (non-hermitian) adiabatic invariants are related to the usual creation 
and destruction operators through a Bogolubov transformation and therefore do not cor¬ 
respond to the latter. In analogy with the particle destruction and creation operators, 
however, one can construct coherent states of the adiabatic invariants. We have also eval¬ 
uated the expectation values of diverse quantities of physical interest with respect to such 
states and shown, for example, that the expectation value of the scalar field operator is a 
solution to the classical equation of motion. 

In order to obtain some estimate of the corrections to the adiabatic approximation 
we have determined the diverse quantities to lowest order in 6 (the ratio of the Compton 
wavelength of the scalar particle to the Schwarzschild radius of the black hole) and for small 
enough times tj (nonetheless satisfying > S). This has shown that the dominant term 
{0{5/rf‘)) comes from particle creation in the average matter Hamiltonian, whereas the 
corrections due to the fluctuations in the number of produced particles, in the uncertainty 
principle or in the semiclassical approximation are of higher order /rf)). To this 

order of approximation (see Eq. (3.45)) there is no distinction between the different three 
geometries (e = 0, ±1). 

Clearly the production of matter will influence the collapse and one may consider the 
associated back-reaction. As we have mentioned, the leading term is of order 0{5/'rf‘). It 
comes from the expectation value of the matter Hamiltonian with respect to the adiabatic 
invariant coherent states and we have seen that it leads to a slowing down of the collapse, as 
expected. It would be of interest to also consider the effect of matter being radiated away, 
which would need the extension of our results to a Vaidya, rather than a Schwarzschild 
metric. 
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